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Abstract 

Recent STM measurements have observed many inhomogeneous patterns of the 
local density of states on the surface of high-T c cuprates. As a first step to study 
such disordered strong correlated systems, we use the BdG equation for the t-t'-t"- 
J model with an impurity. The impurity is taken into account by a local potential 
or local variation of the hopping/exchange terms. Strong correlation is treated by 
a Gutzwiller mean-field theory with local Gutzwiller factors and local chemical 
potentials. It turned out that the potential impurity scattering is greatly suppressed, 
while the local variation of hoppings/exchanges is enhanced. 
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1 Background 



Anderson's theorem tells us that the 
s-wave superconductivity is insensi- 
tive to small potential scattering. On 
the other hand, the d-wave supercon- 
ductivity has zero superconducting 
gap in the nodal direction, and thus 
may be sensitive to disorder. How- 
ever, experimental observation of 
the high-temperature superconduc- 
tivity, which many people are nowa- 
days convinced has d-wave symme- 



try, seems robust against disorder. 
For example, local den sity of states 
measured by the STM (IKohsaka 071 ) 
show clear V-shape at low energy 
that indicates the d-wave nodes are 
not much influenced by disorder. 
Theoretically, it is proposed that 
this protection of V-shape is due to 
strong Co ulomb repuls ion between 
electrons ( Anderson 00l ). Hence de- 
tailed studies of effects of strong cor- 
relation for impurity scattering are 
necessary, and in this paper we focus 
on it. 
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In the Gutzwiller approximation, the 
model parameters are replaced by 
renormalized ones in return for tak- 
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ing the intractable projection opera- 
tor away. Then, i-term is renormal- 
ized by a factor g t < 1 because hop- 
ping is more difficult in the presence 
of projection. On the other hand, 
J-term is renormalized by g s > 1 be- 
cause each site is more often singly 
occupied. In this paper, we focus 
on another term; how are impurity 
terms renormalized? 



(ii) local t variation, 
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(iii) local J variation, 

H imp = SJ So" Sj. (6) 
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2 Model 



3 Method 



We use t-t'-t"-J model with an impu- 
rity term, namely, 



H = H t + Hj + H imp , 



(1) 



H t = P G \-J2tijcUja}P G , (2) 
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where tij = t,t',t", for nearest, sec- 
ond, third neighbors, respectively, 
and otherwise zero. The summa- 
tion in the J term is taken over 
every nearest-neighbor pair. The 
Gutzwiller projection operator Pq 
prohibits electron double occupancy 
at every site. Throughout this paper 
we take t' = -0M,t" = 0.2t, J = 
0.28i, and the hole density x = 0.125. 

We put an impurity at i = 0. Here, 
we try three different types of the im- 
purity term and compare them: (i) 
impurity potential, 



H imp = -Pg (v J2 C la c 0aj Pg, 



(4) 



We solve a Bogoljubov-de Gennes 
(BdG) equation using the Gutzwilller 
approximation with local Gutzwiller 
factors and local chemic al potentials 
(IQH Wang 06l : IC Li od ). Let us as- 
sume that a good variational ground 
state can be represented in the form 
of P\ip), where represents a state 
obtained by solving a BdG equation. 
The projection P includes Pq and 
an operator to control the particle 
number. 



The Gutzwiller approximation as- 
suming that the average of the local 
electron density n. 
changed by the projection, i.e 
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gives the local Gutzwiller factors as 
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where the local hole density Xi = 
1 — (ni)o- Here, (. . .) denotes the 
expectation value by From the 
extremum condition the free energy, 
the following Bogoljubov-de Gennes 
(BdG) equation is derived: 



H 
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where Xij = (cl T c jT ) = (c^c^) , 
Ajj = (cjjCj|) . We assume that x^- 
and Aij are real numbers and that 
Ajj = Ayj. The local chemical poten- 
tial /ij is the derivative of the internal 
energy with respect to the local hole 
density, 
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We use a supercell composed of 
20x20 sites whose origin is an im- 
purity site. This supercell is re- 
peated to construct a superlattice of 
10x10 super cells wit h the periodi c 
boundary condition (ITsuchiura 00l ). 
Then, the Hamiltonian can be block- 
diagonalized by the Fourier trans- 
form with respect to the supercell in- 
dices, and calculation of expectation 
values is reduced to an average over 
many twisted boundary conditions 
of the 20x20 site system. 



4 Results and Discussion 



Let us start from impurity type (i), 
namely, the impurity potential. As 
mentioned in Eq. (7), is not renor- 
malized by any g factor by definition. 
However, as one can see in the BdG 
Hamiltonian, Eq. (10), the impurity 
potential can be compensated by 
Therefore, we define a renormalized 
impurity potential by including dif- 
ference of fj,i, namely, 



V = V - (//»= 
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where /i^ is the local chemical poten- 
tial at the site farthest from i — 0. 

Figure 1 shows Vq function of Vo. 
Note that Vo is strongly suppressed, 
and the renormalization factor seems 
about the order of g t = 2x/(l + x) 
which is plotted by a dotted line for 
comparison. In our understanding, it 
can be explained as follows: Basically, 
impurity sites are uncomfortable for 
electrons to stay at. However, under 
strong electron correlation, every- 
where is uncomfortable and impurity 
sites may be less uncomfortable com- 
pared to the systems without corre- 
lation. This effect appears as com- 
pensation of the impurity potential 
by local chemical potentials. 

Next, in order to see the influence 
from the type-(ii) and (iii) impu- 
rity, we plot the local electron den- 
sity along x-axis in Fig. 2. Here, 
Fig. 2(a) is for the type-(ii) impu- 
rity with 5t/t = ±0.1. When t is 
smaller (larger) locally, the electron 
density near the impurity site is 
higher (lower). Therefore, effective 
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Fig. 1. The renormalized impurity po- 
tential as a function of the bare impu- 
rity potential. 

tij, namely, ghUj, is further smaller 
(larger). This should be because the 
system try to gain energy by increas- 
ing local Gutzwiller factor at large 
t region. Namely, this defect is en- 
hanced by the strong correlation. 

We also find similar behavior for 
type-(iii) impurity as shown in 
Fig. 2(b), where 5 J/ J = ±0.1. When 
J is smaller (larger) locally, the elec- 
tron density near the impurity site 
is lower (higher), and gfjJ, is further 
smaller (larger). However, the mag- 
nitude of the enhancement is much 
smaller than that in Fig. 2(a). This 
should be because dgfj/dxi is much 
smaller than d<7*-/dxj, and the effect 
by local electron density modifica- 
tion is also much smaller in the J 
variation than in the t variation. 
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Fig. 2. The local electron density along 
x-axis when an impurity is located at 
0: (a) a type-(ii) impurity (t- variation), 
(b) a type-(iii) impurity (J- variation). 
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